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^ ■ Abstract 

We study small x divergences in boost invariant similarity renormalization 
group approach to light-front QCD in a heavy quark-antiquark state. With 
CNj ' the boost invariance maintained, the infrared divergences do not cancel out in 

the physical states, contrary to previous studies where boost invariance was 
00 ■ violated by a choice of a renormalization scale. This may be an indication that 

the zero mode, or nontrivial light-cone vacuum structure, might be important 
for recovering full Lorentz invariance. 

^ ■ I. INTRODUCTION 

Qh! Light-cone or light-front description of field theories has long attracted attention for 

I various reasons, often related to its peculiar kinematics 0. Light-cone coordinates are used 
in vastly different applications usually in Hamiltonian formulation (see for an extensive 
list of references), ranging from deep inelastic scattering and phenomenology to attempts to 
formulate the theory of everything. Our own interest in light-front field theories is motivated 
5^ ! by an assumption that in light-front coordinates it is possible to derive from first principles 
a self-consistent and systematic constituent approximation to QCD We will give more 
details below. 

As a consequence of the light-front kinematics, one component of the three-momentum, 
so called longitudinal momentum k^, is positively definite and can be interpreted as playing 
role of a Newtonian mass . Thorn has long advocated that for this reason, light-front for- 
mulation of string theory is one of the best hopes for a truly fundamental, non-perturbative 
description of strings |0 and to quantitatively realize the conjectured correspondence be- 
tween string and field theory 

Because of the above mentioned properties of the longitudinal momentum k'^, light-front 
vacuum was thought to be trivial. This notion turned out to be naive, but the light-front 
vacuum indeed can contain only particles with k~^ = known as the zero mode. Therefore, 
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the light-front vacuum can be made trivial by imposing appropriate cutoffs. This is the 
last point that the light-front community agrees on. Opinions differ when it comes to the 
question how to proceed further. 

Basically, there are two approaches to the problem of light-front vacuum ^j. One is 
to address it head-on 0. Typically, this is done for various, often lower dimensional, field 
theories (for review see |jlO[, for extensive list of references see also 0) in the context of 
discreet light cone quantization (DLCQ) (for a recent review of the method as well as an 
expansive list of references see [|lll). The theory is formulated in a box of a finite longitudinal 



size which makes discreet and non-zero. The zero mode |]T^ is then determined as a 
solution to constraint equations. Constraint and dynamical zero modes are often discussed. 
Critics complain that so far, it is not clear how much of the zero mode determined in this 
manner is an artifact of the method (in other words, a counterterm specific to the DLCQ) 
and whether the continuum limit existsQ. 

As an untypical example of direct studies of the zero mode, we would like to mention 



a recent work by Tomaras, Tsamis and Woodard [Tj] on back reaction in light-cone QED. 
Though motivated by the back-reaction in quantum gravity occurring on an inflating back- 
ground, their work addresses some issues of the light-front vacuum without having to evoke 
DLCQ. They had constructed a full operator solution to a free QED coupled to a constant 
external electric field in continuum (3+1) dimensions. In this set up, all modes are forced 
to go through the zero mode at which point particle pairs are created. The zero mode of 
the constraint components of the fermionic field is shown to be crucial for unitarity. 

The other approach to the problem of zero mode is more pragmatic. Instead of trying to 
solve for the zero mode, it is simply cutoff, be it with DLCQ |T5[ or an explicit infrared cutoff 



in a continuous formulation 0. Physics associated with this mode can then be put in form 
of counterterms, if needed, for example, to restore symmetries or account for phenomena 
associated with the vacuum. Traditionally, spontaneous symmetry breaking was viewed 
as an example of such a phenomenon |]^|T^JT^. However, Rozowsky and Thorn have 
argued recently that, while conceding that the inclusion of a fundamental zero mode is 
a valid theoretical option, it is not necessary to describe spontaneous symmetry breaking 
where its presence seems to be most needed. Indeed, in scalar quantum field theory in 
(1+1) dimensions DLCQ the physics of spontaneous symmetry breaking is completely and 



accurately described without the zero modes |]T5 



The need to put in new counterterms associated with the infrared (IR) regulator in 
our continuum formulation |^ was anticipated |T^ but was not encountered yet in the 



applications to hadronic physics so far |T6|. Perry [|T7| has shown that even though the one 
body and two-body effective operators are each separately divergent as goes to zero, 
the divergences exactly cancel in any color singlet state. The cancelation does not occur 
for non-singlet states, leaving them with an infinite mass. This together with a naturally 
generated confining potential (imprecisely referred to as "logarithmic" ) is a plausible feature 
of the approach. Note that both the effective confining potential and the infinite mass of 
the color non-singlets originate from small regions. 



^See for example, transcript of discussion sections [I3| at the Fourth Workshop on Light-Front 
Quantization and Non-Perturbative Physics. 
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A shortcoming of the above described result is that it was obtained in a similarity renor- 
malization group limited to matrix elements that required a specific infrared regulator (theta 
function) and an introduction of an arbitrary scale which violates an explicit kinematic 
symmetry of light cone. In a bound state calculation one can argue that there is a preferred 
scale, i.e. one associated with the typical longitudinal momentum of the state, or, the total 
center of mass however, consequences of the violation of the kinematic symmetry are 
not known. 

Since then, the similarity renormalization group approach has advanced so that it is no 
longer necessary to violate the kinematic boost invariance, and to generate counterterms 
dependent on the total center of mass . It is also possible to use an arbitrary form of the 



small X regulator [|T^] (x being a dimensionless, boost invariant fraction of the total P^)- 
Thus, for the first time we are able to study with some degree of generality the issues related 
to the light-cone zero mode. Recently, Glazek has found |jT9| that even though infrared 
divergent terms cancel out in the running coupling, there is a residual finite dependence 
on the functional form of the infrared regulator. In this paper, we wish to study the issue 
of small X (or infrared) divergences in color singlet states consisting of a heavy quark and 
antiquark of the same fiavor, for simplicity. 

The paper is organized as follows. Section II is devoted to an overview of the similarity 
renormalization group approach light-front QCD. We start with a general description of 
the approach. For the sake of making this paper self-contained, the reader is reminded of 
the structure of the classical light-front Hamiltonian. The corresponding quantum Hamilto- 
nian cannot be defined without regularization. Details of regularization determine structure 
and form of counterterms. In this case, details of regularization are particularly impor- 
tant, because only one of the regulators can be properly removed by renormalization. After 
introducing regularization, we present a technical description of the boost invariant simi- 
larity renormalization group for particles. Remaining sections deal specifically with quark- 
antiquark color singlets to second order in g. Section III contains the effective one body 
operators to second order in g\ section IV the effective two body operators to the same order. 
In section V we address the issue of infrared divergences in those operators. We conclude 
with a short summary and conclusions. 



II. SIMILARITY RENORMALIZATION GROUP APPROACH TO QCD 

In this section we briefiy review the similarity renormalization group (RG) approach 
to light-front (LF) QCD, introduced in ref. 0. We show the unregulated canonical light- 
front Hamiltonian. Then the regularization that we use is briefiy explained, and finally, the 
similarity renormalization procedure is outlined. 

The basic assumption behind the approach is that it is possible to derive a constituent 
picture of hadrons from QCD. To separate vacuum fiuctuations, it is convenient to use 
light-cone coordinates with cutoffs preventing zero longitudinal momentum. Then vacuum 
is forced to be trivial. In addition, since such a cutoff introduces a nonzero minimum for 
the (kinematic, positive) longitudinal momentum, it also restricts number of particles in 
any state with a fixed longitudinal momentum P^ . For massive particles, due to light-cone 
free energy increasing with decreasing longitudinal momentum, many-body states tend to 
have higher free energy than few body states. Mixing of low energy few particle states with 
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many-body states can be expected to be naturally small, at least at weak coupling. These 
features make for useful prerequisites toward the constituent picture of hadrons. 

The apparent difficulty with renormalization of light-front Hamiltonians (compared to 
Lagrangians) is turned into an advantage by using similarity renormalization ||20|. The basic 
idea of the similarity renormalization group is simple. The regulated bare Hamiltonian that 
mixes all energy scales is transformed via a unitary transformation to a Hamiltonian that 
contains direct couplings only between neighboring scales. Such a Hamiltonian is referred to 
as "band-diagonal" in the sense that it does not allow for direct couplings between arbitrary 
scales. At any finite order of perturbation theory, a band-diagonal Hamiltonian cannot 
produce any ultraviolet divergences providing its matrix elements are finite. Therefore, by 
requiring that the band-diagonal Hamiltonian be independent of the regulator, counterterms 
that need to be added to the bare Hamiltonian can be identified. 

Upon completion of renormalization, the effective band diagonal Hamiltonian is subject 
to diagonalization. We refer the reader to [TE\ for more details regarding this step of the 
procedure. 



A. Classical light-front QCD Hamiltonian 

The starting point is the canonical light-front QCD Hamiltonian in light-cone gauge, 
A'^ = 0. We will not explicitly show terms that are not important for the specific calculations 
presented in the next sections. For a detailed discussion of the light-front Hamiltonian see 



Purely gluonic terms are studied in [|T^] 



The part of the classical unregulated canonical Hamiltonian that is relevant for our study 
of the bound state of a quark and antiquark of the same fiavor is, 

H = Hf,,, + V, + V2 ; (1) 

where iffree is the free light-front Hamiltonian: 

i//,ee = y dx-d\^ f -7/;7+ ^ - -^[d^A- - d-A^)[d^A, - d,A^)\ , (2) 

and 

Vi = g dx~ d^x±il>^il) (3) 



contains the standard order g quark-gluon coupling. Here if) and A^ = X)a ^a^" ^ire free 
light-front fields: 



and 



= ( = 0, A- = ■ A^, ) . (5) 
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The constrained fields contain, in addition to if)- and A~ shown above, 0{g) terms that are 
fully determined functions of the physical degrees of freedom. These produce new terms in 
the canonical Hamiltonian, among which 

V2 = -2g'J da;-rfV(^lT>+) (^iT'^^+) (6) 

is the so-called instantaneous gluon exchange between two fermions. As the next step, fields 
are expanded in mode functions consisting of (light-cone) spinors Upa,Vpa {p denotes 
momentum, a denotes spin), SU(3) color spinors Xc and plane waves, 

^ (x", x^) = E / {2^f2p+ [^-"P-^"'"" ^P'-.- + XcVpae'^^'dl,^}^ . (7) 

The Hamiltonian is expressed in terms of what will become, after quantization, creation and 
annihilation operators of fermions {h\ d) , b, d) and gluons a^, a. For example, the kinetic 
part Hfree bccomes 

(27r)32p+ p+ W-^P- + 4cAc<x)+Ey ^^^^)i^-^<<x«P- • (8) 

where Ep = ^ ^ and Ek = are free light-front energies for massive and massless 
particle, respectively. 

As another example, consider the interaction V2 that contains (among other terms). 
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{2ixf2p^ 

~^PlCio-i'^P3C30-3^P2C20-2'^P4C4CT4 + ••• (9) 



where Xi = pf/P^ is the fraction of the total carried by the particle i. Momentum 
conservation is assumed implicitly. Similarly, Vi is also rewritten. 



B. Regularization 



As will be illustrated in Section III, the ultraviolet (UV) and IR divergences in the light- 
cone QCD are interwined (for additional discussion see Apendix E of ref. |19[)- Fot this 
reason and because at present only UV renormalization is achieved in a systematic manner, 
details of regularization are important. We regulate the canonical quantum Hamiltonian 
with cutoffs on the changes in transverse momenta and on longitudinal momentum fractions 



as in ref. |jT9[, term by term in expansion of creation/annihilation operators. For example, an 
operator consisting of one creation and two annihilation operators would contain product of 
two regulators, one for each of the annihilation operators. Regulation of operators containing 
more than one creation operator (and their Hermitian conjugates) is slightly more involved. 
For details see |]T^. The instantaneous interaction (^ is viewed for regulation purposes as 



an exchange of a virtual particle between two vertices. 
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These general rules can be satisfied by various cutoffs. Specifically, we use a regulating 
factor of the form (arguments are generic for the moment) 



(k^^^ x) = exp(-fc^ VA')r5(x) . (10) 



that is particularly convenient for analytical calculations. Here and in what follows is an 
infrared regulator that prevents its argument to become zero. We wish to study dependence 
on this regulator, hence we leave its functional form unspecified. A is an ultraviolet cutoff. 
The limit with the cutoffs removed is achieved by A ^ oo, 5 ^ 0, in which case r^^ 1. 
In order to preserve kinematic symmetries, we choose the arguments of the regulators to 
be kinematically invariant quantities, in particular, Jacobi momenta defined with respect to 
the vertex in question. As an example of a regulated interaction, the standard quark-gluon 
coupling has the following form in the regulated quantum theory. 







12 X2\ 
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3^3 / 



U2 {eTl>.) U3T^%^ bla\h + h.c. + ...j . (11) 

For clarity, summation indices are implicit here, sum is over spins and colors of all 3 particles, 
integration runs over all 3-momenta with a momentum conserving S function, and index i is 
short for {pi{ki), (Xj, Cj) as appropriate, u, u are light-cone spinors, is a gluon polarization 
vector, = (e^ = 0,e'j^^ = 2k^-e^/k~^,e^], and here al = al^o-ici denotes specifically 



a gluon creation operator. (We did not intend to list all terms in Vi). The longitudinal 
fractions are Xi/x^ = k^/p'^, Xi/x^ = kf/p^, satisfying 

Xi + X2= Xs, (12) 

and the transverse Jacobi momentum Kj;2 is defined by 

fci = 59o -|- Kin 

x/^ 

P2 = —pi - ■ (13) 

With the specific choice of the functional form of the UV regulator in (0), the entire 
regulating factor in (|1T[; 



IS 




exp|-2^1..(|M|). (14) 



This concludes our definition of the bare regulated Hamiltonian. 



C. Boost invariant similarity renormalization 

The dependence of the Hamiltonian on the ultraviolet (UV) regulator A is removed via 
similarity renormalization scheme whose basic idea has been outlined in the introductory 
part of this section. Similarity renormalization can be performed in various ways, in terms 
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of matrix elements (the original Glazek- Wilson formulation) or in terms of coefficients of 
creations and annihilations operators (Glazek's similarity for particles), and for various 
definitions of the similarity formfactors in the effective Hamiltonian [^8|,^. For a review and 
comparison of numerical efficiency of various formulations for Hamiltonian in QCD in (2+1) 
dimensions see p3[. The focus in this paper is restricted to the similarity renormalization 



group for particles (i.e. for creation and annihilation operators) that we use [T^. For 
introduction to the similarity renormalization group for particles see the latest reference in 

Let us denote the regulated canonical bare light-front Hamiltonian, -f^canAS, together 
with the as of yet undetermined counterterms, Xa, in accordance with |]19[ as H^s- H/^s is 
written in terms of creation/annihilation operators of bare, or current, quarks and gluons. 
For the purpose of this discussion, we denote all creation/annihilation operators, bosonic 
and fermionic, generically a\ a. If this current Hamiltonian were to be used to describe 
mesons, physical states would have to be very complicated superpositions of current quarks 
and gluons. We wish to find an effective Hamiltonian that provides a simpler description at 
hadronic scales A. To arrive to a simpler picture of mesons, the effective Hamiltonian should 
rather be expressed in terms of constituent quarks and gluons. Yet, the two Hamiltonians 
are equal, 

H^sia) = Hxiax), (15) 

meaning that they have the same eigenvalues. A is a parameter of dimension of mass that 
will be explained in more detail below. As the Hamiltonians are related by the unitary 
transformation Ux, so are the current and constituent particles, 

ax = UxaUi. (16) 

The creation/ annihilation operators , a satisfy the usual fermionic or bosonic commutation 
relations. We use Lorentz invariant normalization 2(27r)^p"^5^(p' — p) for the commutation 
relations. 

Applying the transformation Ux to both sides of (|15D and owing to its unitarity, one 
obtains 

UlH^s (a) Ux = Hx (ulaxUx) (17) 



From (|T^ follows that all dependence on A in ( p!7| ) is in the coefficients of the various 
combinations of creation/annihilation operators. So differentiating with respect to A affects 
only the coefficients, not the operators, and we obtain 

:^^a = -[t-a,Ha] (18) 



dX 



where we have denoted 



Hx = Hx [UlaxUx) = Hx (a) (19) 
in accordance with [|T9[, and tx = UlU'x generates infinitesimal transformations. Equation 



(|T8|) is an operator equation, but apart from c-numbers it contains only known, current. 
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creation/annihilation operators and their commutators. For this reason, it can be solved in 
an operator form, in contrast to earher formulations in terms of matrix elements. 

The transformation is constructed from the requirement that the Hamiltonian H\ be 
band diagonal, which means, as we stated in the introductory part of this section, that it 
does not contain direct couplings between arbitrary scales. Parameter A is a measure of 
width of momentum space formfactors fx that appear in vertices of the renormalized H\. 
For this reason, A is often referred to as band width. Schematically, we require that 

Ha = fxGx (20) 

which indicates that we want Tix to vanish in the part of the phase space were the similarity 
formfactor fx is zero. This equation also defines GaQ- 

It is easier to work with Gx than with T^aQ- Gx is split into the free part Gqa = Gq 
that is, in present work, independent of interactions and does not change with A, and the 
interacting part Gjx- Rewriting the differential eqn. ([18|) gives 



^Gjx + h-^Gjx = -[rx, Go] - [tx , Gjx] (21) 

The left side of the equation is independent of Gq because Gq does not change with A and 
because fx does not change on the diagonal. Gjx consists of the canonical interactions as 
well as new effective interactions ensuring that the effective Hamiltonian is equivalent to the 



original current one. Therefore, (21) describes two unknowns, tx and Gjx- We have freedom 
to arbitrarily split this equation into two. Without loss of generality, we can assume that 
when the interactions vanish, so does the generator tx- If we choose 

/a^G/a = - fx [rx , Gix] , (22) 

then 

[rA,Go] = ^((l-/A)G,A), (23) 

and the generator of the transformation is of order Gix- Changes in Gjx, corresponding to 
new effective terms, are of second order in interactions. 

Equation ( ^31) can be solved for tx, for example order by order in the interactions. In fact, 
with Go being the free Hamiltonian (||), the solution is simple. It is the same operator as the 
right hand side of the equation (|23|) with an additional factor [J2a Ei — I]c -^j] ^ ■, where a, c 
stands for annihilated/created, and i?'s are free light-front energies created/annihilated in 
the vertex (see (H)). The additional factor arises from the free energy in Gq and contracting 



Be aware that Glazek in |1S] uses notation Qx to make a distinction with his earher works. 

^Factorization of the overall /a, and writing the renormahzation group equations for Gx rather 
than Hx is not necessary, but it sure makes the equations much easier to work with. For comparison, 
see the earliest reference in [Efll. 
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one creation/annihilation operator with the conjugate in Gq. It is, however, not necessary 
to exphcitly find tx- If we denote the solution to ( p3l) 



Go 



then the solution to ( P^D subject to the boundary condition 



G 



A— »oo 



AS, 



(24) 



(25) 



where H^s includes both known canonical terms and unknown counterterms, can be written 

as 



Gy 



H, 



can A(5 



+ Xa+ / da 



Go- 



(26) 



Counterterms should remove any dependence on A arising from the second term. This 
criteria is sufficient to determine X^. Finite parts of the second term constitute new effective 
interactions and/or modifications to canonical terms. Finally, the Hamiltonian is found by 
substitution to (pO|) . 

The equation (^) is an operator equation but the creation/annihilation operators are 
not affected by differentiation. They merely play a role for "counting" purposes. If we are 
interested in writing down a differential equation for a coefficient of a specific combination 
of creation/annihilation operators, to any given order in g we can see what operators in the 
interaction Hamiltonian need to be included on the right-side of the equation, and which 
of the creation/annihilation operators on the right side need to be contracted. Note that 
since the right side of (^) is proportional to a commutator, the coefficients in the effective 
Hamiltonian arise only from connected terms. 

The preceding discussion was valid for an arbitrary choice of the similarity formfactor, 
and we have not specified how the formfactors are imposed. Specifically, we require that 
any combination of creation/annihilation operators in the Hamiltonian is accompanied by a 
formfactor fx; a particularly useful choice of the similarity formfactor fx has the following 
form 



fxiMl~Ml)^fx{ca)=exp 



[Ml-M 



212' 



(27) 



where Ail = createdPif IS the square of sum of all free four-momenta created in the 

vertex and similarly, M.^ = annih.P'j) is the square of all free four-momenta annihilated 
in the vertex. 

For example, for any finite A Hamiltonian, the integrand of the operator ([Tl|) is modified 
by insertion of a factor 



f\{{ki+P2f -pI) = exp 



P1+P2) 



212' 



A^ 



(28) 



The corresponding term in Gx, i.e. the Hamiltonian without the vertex formfactor, is to the 
lowest order in interactions just equal to ( [Til) and independent of A. 
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In this work, we expand the Hamiltonian relevant for a quark-antiquark state in powers 
of g, and find effective interactions up to g"^. The perturbation theory for Gx is well explained 



in section 3.c of reference |T^. As it is often the case, in practice the most convenient way 
to determine the factors arising from the similarity transformation from infinity down to the 
finite scale A in eqn. ( p6D is with the help of diagrammatic rules. Unfortunately, these seem 
to be the best kept secret in the field, and cannot be found listed in any reference, perhaps 
for fear of confusion with the ordinary perturbation theory for the S-matrix |]TB|]. The 
calculation of effective interactions to 0{g^) is much simpler than at higher orders, because 
it involves only G\ at 0{g) which is independent of the similarity scale. Consequently, all 
dependence on the scale in eqn. ( pUf ) is in f„ and df^/da. This is not true in general, 
when care must be taken to keep track of all a dependence. Therefore, the second order 
calculation does not require "full-blown" diagrammatic rules and we will not develop them 
here. We will, however, show the simplified rules sufficient for the present calculation, with 
the understanding that they should and can be generalized for higher order applications. 
To draw a diagram for the effective coefficient of any operator 0\ in Gx-, 

• draw an external line for each creation and each annihilation operator in Ox- Draw 
all connected diagrams connecting the external lines. (At second order, one-body 
operator has only one diagram. Fig. 1, reminiscent of a self-energy diagram in the 
ordinary perturbation theory. The two-body operator has two.) 

• Let us denote total momentum annihilated/created in a vertex i clS Cli , Ci . (At second 
order, i = 1, 2.) Let pt,pf' be the momentum conserved in the vertex i, and light-cone 
free energies = {-Ml. + pi^)/pt : created, and Ea^ = (A^^i + Pi'^)/pt y annihilated. 

At each vertex i, include the standard factors, as well as regulators, already present in 
the Ga, such as, for example, the familiar Uie^'j^u,with a momentum conserving delta 
function, at each qqg vertex i. 

• The novel element comes from the similarity formfactors and scale. The entire diagram 
is to be integrated d a. At second order, assign 



(29) 



to the diagram. This factor contains all the dependence on the scale, and, for f„ given 
in (^), it can simply be integrated, resulting in 

'pt{Ml-MlJ+pt{Ml-Ml 



{Ml - Mir + {Ml - Ml 



[fx{ciai)fx{c2a2)-l]- (30) 



Diagrams can acquire an overall sign and/or combinatorial factors so that the same 
ordering of uncontracted creation/annihilation operators on both sides of the eqn. (^61) 
is achieved HTHl- At second order, however, these overall factors are obvious. 



• Generalization of the procedure so that it is applicable for higher order calculations is 
st raight forward . 

As a specific example, we present in the next section the renormalization of the one body 
operator. 
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FIGURES 




ai 



Ci 



C2 



FIG. 1. Diagram for the quark one body operator at second order in (7, not to be confused 
with the ordinary perturbation theory for S matrix. Thick hues correspond to uncontracted cre- 
ation/annihilation operators. Double hue is a contracted gluon, single thin line is a contracted 
quark. Vertical dashed lines show which momenta add up to Oj, Cj, respectively. All other notation 
is explained in the text. 



III. RENORMALIZATION OF THE ONE-BODY OPERATOR 



The one-body operator was previously derived, but not further studied, in ref. [|TH| with 
slightly different choices. 

The one-body operator when A = 00 is 



d^p 



+ ^ 5m\ ^ 



P 



+ 



(31) 



where is the bare mass, A is the UV regulator and 5m\ is the counterterm to be 
determined so that renormalized Hamiltonian is independent of the regulator. 

When A is finite, there are corrections (with coefficient 6m\/p~^) to the one-body opera- 
tor0 which start at O {g'^)- 



^With boost invariant similarity, corrections have the proper dispersion relation. 
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To visualize the diagrammatic approach outhned at the end of previous section, we 
present in Fig. 1 a diagram for the lowest order correction to the one body operator with 
the notation used in the list of rules. For example, momentum oi is the free momentum 
of the uncontracted quark p; momentum Ci in this diagram is equal to sum of the free 
momentum of the virtual gluon kg and the virtual quark kg. Three momenta are conserved, 
so only one of the virtual particles momenta is independent. 

One factor of ( [TT| ) (with appropriate momenta and without integration over the external 
momentum) stands for each vertex. Sums of the polarisation vector, spin and colors of the 
virtual particles lead to standard expressions involving color factors and matrix elements of 
Dfj,u{kg)'~f'^ {j^q + m) 7*^ in the spinors of the uncontracted quark. It is convenient to express 
the virtual momenta in terms of Jacobi momenta K"*", x. This loop momentum is integrated 
over. 



The similarity factor (|3 



further simplifies. Using that pt = pt = P' 
^/(l — x) in this case, we get 



Ml 



m 



Mi = Mt 



and Ml^ 



rn? — 



^ l-x) 



f. 



2 + 
m — ( 

X 



+ 



1 — x' 



(32) 



In order to write the one body operator with the same ordering as the original mass term, 
the uncontracted 6, W have to be commuted which gives an overall (-1). 

After this straightforward algebra, we find that the correction to the one body operator 
is indeed of the expected form constant / p^ b'^b (independent of p-^) with the coefficient 



exp 



' A2 



2y + 4 



1-y 



Am? 



z + -r—vrfi 



y 2\2 

m 



(33) 



independent of p^. Here m = moo; z arises from integration over transverse momenta and 
carries a dimension of mass squared. The coefficient ( ^31) consists of terms /i^ independent 
of A, containing UV/IR divergences, and A-dependent terms, both finite /i| and IR divergent 



(47r)2 2N 



2 I 2 I 2 



(34) 



where: 



= -A^ cs + 4m2 log — - 4m2(7 + log 2 - 2) 



- A^v^ 



y 2\2 



\^V2^ dy 
Jo 



l+/'d,i^Erffv/2^^-i^ 

Jo y V ^ 1 - y 

rsiy) 




y 



4m 



z + -^m? 



(35) 



(36) 
(37) 
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Note that the coefficient of the quadratic divergence = Jq dy[r1 
diverges as 5 ^ 0. 

The one-body operator can be written as 



j{l-y)/{l-y) + 2r]{y)/y'] 




(38) 



where 



(47r)2 2N ^ 



(47r)2 2N ^^'^ 



(39) 



The UV divergent fij^ is combined with the counterterm 6mj^ so that the dependence on 
the regulator A is removed. Nevertheless, the mass m\ is not finite in the limit 6 —>■ 0. In 
what follows, we are concerned only with the behavior in this limit, rather than details of 
the finite terms. For this reason and for simplicity, we grouped all finite terms contributing 



We are seeking a coefficient V of an operator consisting of two creation and two annihi- 
lation operators. 



which we expressed in terms of relative (Jacobi) momenta and set the total transverse 
momentum of the initial and final two-body states to zero because of kinematical boost 
invariance. We also factored out the total because all regulators and form factors are 
expressed in terms of boost invariant quantities. The conservation of the three momentum 
is implicit here and in what follows. 

The coefficient Va depends on momenta, and contains terms of 0{g'^) and higher We are 
interested only in the color singlet and most IR divergent part of this operator. (This part 



We factored out the omnipresent A^JliiX2X^xl because it cancels with the similar factor in the 
definition of the bound state wavefunction. 



to m\ into m^(A). 



IV. TWO-BODY OPERATORS 




(40) 



that creates/annihilates quark-antiquark pair with free momenta 




(41) 
(42) 
(43) 
(44) 
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is diagonal in spins). Depending on its degree of divergence, it can be concluded whether or 
not the subleading terms are IR convergent. 

The canonical Hamiltonian already contains a two-body operator of this form, the in- 
stantaneous gluon exchange V2 ® . With the infrared regularization as in |jl9[ its coefficient 
is 



e{xi - X3) rs 



X3 

Xi 



(xi - X-iY 
Xi — X3 



+ Oi^X'i - xi)rs ( — ) ^"5 



Xi 



Xi 
X-i - Xi 

X3 



X2 

X4 



rs 



X4 



f xi - x^ 

\ X4 

X2 



(45) 



where rs is the (unspecified) infrared regulator. 

With the choice of fx ( pTj) and regularization as in ref. [|19|, the spin independent, most 
infrared divergent, color singlet part of the coefficient of the effective two-body operator to 
lowest (i.e. second) order in coupling is 



g^Cp (Vi + V2) 



(46) 



where 



6{xi - X3) 



X1-X3 

Xl - X3 



rs(— ] rs 

\Xi/ \ Xl 
-g^CFV2 = -g^CF 



X2 



Xl - X3) 

Xi{Ml^ - Ml) + xi{Ml^ - Ml) 
[Ml,-Ml)^ + {Ml,-M\)^ 

Xl - X3' 



X4 

3^3 - Xl) 
{Xi - 0:3)2 1 (X3 - Xl) 



MMl,-Ml)MMl-Ml)-l 

X3 - Xl 



X2{M% - M 



x3{Mj,-Ml) 



{Ml - Ml)^ 



(Ml, - M^^ 



X3 



X4 



X3 - Xl 

X2 



(47) 



(48) 



It arises from contracting one pair of creation and annihilation operators of a gluon with 



^5 — (^5 /X5, X5, — K,i2 



(49) 



where X5 = 
defined as 



\X3 



Xi\ corresponding to the two time-orderings. Invariant masses M are 



Ml = (pO^ = m 



(50) 
(51) 



at each vertex and rs are the infrared regulators, in accordance with |jT9|. Vi corresponds 
to particle 1 emitting the gluon that is absorbed by particle 2; in V2 particle 2 emits gluon. 
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particle 1 absorbs it. The second line in both expressions comes from integrating the width 
parameter in the similarity RG procedure from infinity down to its value A. 
In terms of the Jacobi momenta, the arguments of the formfactors are 

-2ni,.ki + ^{m' + Ki^) (52) 

+ 2^/€^2 ■ ki + -(m^ + K^^) (53) 

-2K^,.ki+'^{m' + K^^) (54) 

+ 2^/€^2 ■ ki + ^(m^ + Ki^) (55) 



Ml, 


-Ml = 


k\^ 

X2 — 


Ml, 


-Ml = 


Xl 1^, 


X3 X, 


Ml, 


-Ml = 


kP 
xi — 

X5 


Ml, 


-M\ = 


ry.2 1,1.2 

X2 "^5 


Xi X, 



V. INFRARED DIVERGENCES 



The coefficients of the effective one-body and two-body operators diverge when the in- 
frared regulator is removed, i.e. rs —>■ 1. In previous calculations [17| the infrared divergences 
were found to cancel in the color singlet qq states, leaving however, the nonsinglet states 
with infinite mass. The remaining part of the two-body potential was found to be confining. 
A similar observation was made recently regarding QCD in (2+1) dimensions [^]. We want 
to check whether the cancelation occurs here. 

It is easier to address the infrared divergence issue in the framework of a bound state 
equation, or equivalently, in terms of expectation values of the effective Hamiltonian. In 
that case, the creation and annihilation operators are contracted, and the integrands of the 
momentum integrals are just c-numbers consisting of the coefficients of the operators and 
wavefunctions of the bound state. The bound state eqn. reads 



4m^ + AmE]^ 



'12 



(47r)2 N 



1 1 ^^ 

— + — $12 

Xi X2, 



47]-2 



dx'^d^/t 



34 



IT 



[Vi + V2 + V,^,,]hiMl,-Ml,)^ 



34 



(56) 



where (4m^ + AmE) is the eigenvalu^, and = ^(xj, k^) is the bound state wave function 
of a quark- ant iquark state of total momentum P = (-P^? -P"*" = 0^: 



1^) 



dxid'^K^-j 
2{27r)^XiXj 



,/x-x-<^{xi,Kj-.)bld]\0) 



(57) 



^We use this definition of eigenvalue for its convenience in nonrelativistic limit, when 4m^ cancels 
with the same expression on the other side of the equation, and after dividing by 4m the bound 
state equation reduces to the usual Schrodinger equation for eigenvalue E. 
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The first line of tlie riglit side of the bound state eqn. ([56| ) comes from one-body operators, 
the integral part arises from the two-body operator. 

The infrared structure of the expectation value of the two-body operator depends on 
the similarity formfactors in the various coefficients V. The coefficients also contain two 
different sets of infrared regulators corresponding to the two time orderings. To make the 
dependence on all formfactors explicit, let us introduce vi, V2 as 



Vi= fx{Ml,-M\)fx{M 



2 
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f,{Ml-Ml)f,{Ml-M\ 



9{xi - xs) vi 

6{x-i - Xi) V2 



(58) 
(59) 



Vinst does not contain any formfactor, but it does contain the two sets of infrared regulators 
for xs < Xi and x^ > Xi. It is convienient to introduce finsti, 'i^inst2 in analogy with vi, f 2 : 



inst 



e{xi - X3) Vir^tl + 6{X3 - Xi) t;inst2 



(60) 



This defines f's. Note that the infrared regulators in vi and finsti are the same; so are in V2 
and Vinst 2- 

To extract the infrared structure of the integral in (|56|), we first note that the IR di- 
vergence in the coefficients (^51) and (|47|) , (^81) occurs when X5 — and that the small x^ 
behavior of Winsti has opposite sign compared to Vi. For this reason we first regroup the 
integrand in (^) by adding and subtracting finsti multiplied by the same formfactors fx as 
those multiplying vi in (p3), (p3). This splits the integral into two parts. 



471-3 



Cf J dx3d^4^ [Vi + V2 + Vinst] /a(-M?2 - Ml^)<^u = 1 + 1 



(61) 



where 



[e{x^ - X,) {v^ - i;i„sti) [/a(A<35 - Ml) h{Ml, - Ml) - 

+e{x, - Xi) {V2 - t;inst2) [/a(-M?5 " K) fx{M% - Mj) - l] } 

fxiMl2-Ml,)^s, , 



(62) 



and 



47^3 



Cf J dx^d'ki {9{x, - xs) VinstifxiMl, - Ml) fx{Ml, - Ml) 

+e{x3 - xi) Vinst2fxiMl, - Ml) fxiMl, - Ml)} 

fx{Ml,-Ml)^u 



(63) 



Assuming that the wave function $34 is bounded, it is straightforward to show that / 
(|6^) (schematically, the integral containing (1 — /^)(Vinst — Vi — V2)) is bounded and thus, 
has to be infrared convergent. 

The remaining integral, I (|63|) is divergent as X5 0. This means that X3 xi and 
X4 — X2. Upon inspection of the integrand one can see that the formfactors in (^) restrict 
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more and more severely as 0. It appears that, consequently, $34 $12 as the 

divergence X5 = is approached. So we add and subtract $12 to $34, i.e. replace $34 in (|63D 
by ([$34 — '^'12] + ^12)- Then by expanding $34 in Fourier series around $12 it can be shown 
that the part of the integral with [$34 — $12] is, indeed, convergent as expected. Therefore, 
the infrared divergence is contained in 
-,2 



X 



9 



Cf / dxgd^^^ ^(xi - X3) t;i/A(A^^5 - Ml) UMl, - M 



+e{x, - X,) vMMl - Ml) fx{Ml, - Ml 

h{Ml, - Ml,) 



(64) 



multiplying $12. This is similar to the divergent term from the one-body operator. 

To compare the divergence in the two-body operator with that of the one-body operator, 
we need to isolate it. It is sufficient to consider just series expansion for X5 ^ 0. Keeping 
leading and subleading terms, arguments of the formfactors (^) reduce to 



Ml,-M 



Ml 



Ml 



X2 - 2X5 



''12 '^5 



+ 0{xl) 



Ml 



Ml - Mi 



X5 

X5 



Xi 



''12 '^5 

h.1.2 

K,5 



+ 0{xl 



X. 

X4 



2 _|_ 2*^2 1^12 ■ ^5 



X4 



The overall similarity formfactor is in this limit 



h.1.2 



+ 0{xl) 



[1 + 0{x,)] 



(65) 
(66) 
(67) 
(68) 

(69) 



h{Ml-Ml) 

The formfactors contain factors k^"^ /x^ which make the limit X5 somewhat obscure. 
This can be evaded by a change of variables k^'^/x^ = /C^^. In these variables, 

da;5 



X oc 



X5 



[6{xi - X3)rs... + 6{x3 - xi)rs...] d K, 



[1 + o{^5)] 



(70) 



The subleading correction gives / (ix5{y/xl) which is finite. The integral over /C can be 
done analytically. Out of all infrared regulators, only those regulating X5 — > are necessary, 
because the integral only diverges at its lower limit. The remaining infrared regulators can 
be set to 1. The result is 



9 ^ 



I ^2 



' —rs {xi y) rs (xa y) 
y 



—rs (xi y) rs {x2 y) 

y 



or, 



4^' 2 ^a;2^a; 



2 / —rs (xi y) rs {x2 y) + log f^— 
y \X1X2 



(71) 



(72) 



If this term multiplied by $12 is subtracted from the integral in (|56D, it becomes finite in 
the limit 5 — > 0. 
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A. Cancelation of IR divergences? 



Let us summarize the infrared divergent terms in the bound state equation 
We have rewritten the bound state eqn. (|56D as follows 



12 



+ W?{\) 



47r2 



X1X2 

^2 



12 



dxsd^/t 



[Vi + V2 + Vin.t] /a(A<?2 - -^34)^34 



X2 



V2 X1X2 JO y 



dy 

/ — rs{xiy)rs{x2y)^i2 
Jo y 

idy 



X 



^5 (xi y) rs (x2 



12 



(73) 



If the last line of (|73D vanished, the bound state equation would be independent of 6, just 



as in previous formulation of the similarity renormalization or, as recently reported, in 
(2+1) dimensions p3|. 



It is obvious that this does not happen here, except for the leading order in the non- 
relativistic expansion. Indeed, in the nonrelativistic limit Xi = 1/2, and the two different 
functions of xes, constituting the coefficients of the infrared divergence from the one-body 
and from the two-body operators, have the same value (4/^/2) at this point. However, an 
arbitrarily small deviation from x = 1/2 introduces a positive divergent constant into the 
bound state equation in the limit 6 0. The bound state equation therefore is not defined 
in this limit. 

For any finite 6, the bound state equation is well-defined and can be solved. If 6 is small, 
the IR divergence forces the wavefunction to be peaked at x = 1/2. Detailed behavior 
of the wavefunction is likely to depend on the specific form of rs. In any case, the wave- 
function width determines the allowed range of X5. There is no such direct mechanism to 
restrict the transverse momenta, but they are pushed to small values by the formfactors in 
[Vi + V2 + Vinst] as X5 is pushed to zero by rs- It is peculiar that while the similarity formfac- 
tor f\{J^i2~ -Msa) also restricts the range of transverse momenta, for X5 = this restriction 
is less severe than for any other possible X5. This means that the implicit formfactors are 
more important when X5 ^ 0. 

If we assume that 



Xl 



1/2 + e. 



(74) 



with |e| << 1/2, (consequently, X5 << 1/2) and that other momenta, compared to mass, are 
also small, the bound state equation (^) can be expanded. Care must be taken, however, 
to keep track of the infrared divergence. We find 



AmE $12 



+ Ae^m^ + 



m2(A) 



Ae^m^ 
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47r2 



34 



vr 



1 — exp 



1 



y 



exp 



+ Axlrri^f^ 
2xiX^ 



X5 

2 



47r2 



V2 



Jq y 



12 



(75) 



Here E is the eigenvalue shifted by a constant from the one body operator, E = E — 
(m^ — m?{X)) /m. On the right hand side of the equation, the third term is a correction to 
nonrelativistic kinetic energy. It is hkely that this term is smalL The difference between 
the current mass and mass including finite A-dependent corrections can be expected to be 
negligible for heavy quarks. The last expression in curly parenthesis is simply the subtraction 
of the IR divergence in the two body operator to this order in momentum expansion, so that 
the entire potential is final. The last line contains the residual divergence for x 7^ 1/2 to its 
lowest nontrivial order, O (e^) . 

The structure of the right hand side becomes more clear upon substitution 



2x5m = q^; 
$12 = ^(p); 



k' 



"12 



$34 = ^p + q] 



In these variables, 
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The bound state equation (ITH) reads 

- m^(A) 



2p- q + q^ 



(76) 
(77) 
(78) 



(79) 
(80) 



AmE(^{p) 



3_ 

A-K^ 

1 



p 



dx-jd^K 



-Pz 



34 



vr 



Am' MMl,-Ml,)^p + q) 



2^4' 



2m q 

V2 





1 




+ 




exp 


2 




dy 
— r 




Jo 


y 



2^4' 



2m q 



{- 

\Am 



+ 



47r2 



72 



{Pz 



m^ 



y 



rsiy) 



Hp)- 



Hp) 



(81) 



The last line in (|8ll) is the uncanceled divergence which prevents us to take 6 = limit at 
the Hamiltonian level. Without numerically solving (^l]) we cannot make any quantitative 
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statements, but some conclusions can be drawn just from the structure of the bound state 
equation. If not for the Coulomb potential, i.e. the first integral term in (0), it is obvious 
that with decreasing S, the wavefunction would get squeezed more and more. In the limit, 
the wavefunction would basically collapse to the point x = 1/2 and zero relative transverse 
momentum. This might allow the eigenvalue to stay finite, and possibly roughly the same. 
With the Coulomb potential in place there could be some interplay between the positive 
6 divergence and the Coulomb. Note, however, that the Coulomb potential in our bound 
state equation is incomplete due to the presence of the exponential. The exponential cuts 
off the Coulomb potential at small q (and nonzero g^). It is still possible (although not 
guaranteed) that the mass of the eigenstate is finite and convergent in the 5 — > limit, but 
at the expense of an unphysical wavefunction. 

VI. SUMMARY AND CONCLUSIONS 

Owing to nonzero mass of quark/antiquark, the only infrared divergence in the one- and 
two-body operators in a quark plus antiquark state is due to 0, i.e. the virtual gluon 

longitudinal fraction going to zero. 

The infrared divergent part of the one-body operator obeys the one-body dispersion 
relation by construction of the similarity transformation. Specifically, it goes like a con- 
stant multiplying x^^, and is independent of the spectator. The expectation value in a 
quark/antiquark state is a simple sum of the quark and the antiquark contributions. But 
the IR divergence from the two-body operator is a different, more complicated function of 
xes. In addition to an overall explicit function, multiplying the divergent integral, whose 
form can be traced to the choice of fx, there is an implicit dependence through the arguments 
of the regulators. Thus a different x-dependence can be generated for different regulators, 
contrary to the one-body operator. This occurs because boost invariance dictates the form 
of the one-body operator, and at the same time the two-body operator cannot lead to a 
function of the same form. 

In the leading order nonrelativistic reduction, all dependence on infrared regulator drops 
out. Beyond nonrelativistic leading order the divergences do not cancel out in general, 
and it is not obvious that a regulator such that the divergence would cancel out for all x, 
exists. Caution is in place when the nonrelativistic limit is considered, because the infrared 
divergence, even though related to the most infrared divergent part of the effective terms, 
enters at the subleading order in the momentum expansion. 

What are physical implications of this result and how to evade the problem? This is a 
subject of an ongoing study. Here we only outline some possibilities. We wish to emphasize 
that at this point they are just speculations. 

The first question one needs to address is whether the IR divergences are to be cured at 
the level of Hamiltonian. In previous section we argued that it is possible to solve for the 
bound state with any fixed value of S and at the end look how the solution behaves in 5 — >■ 0. 
The eigenvalue can be expected to be finite and convergent, however, the wavefunction in the 
limit must collapse to a point. Such a wavefunction is counterintuitive. It is not necessarily 
a no-go, however. Wavefunction by itself is observable only through expectation values of 
other operators. These operators have to be regularized consistently with the Hamiltonian. 
Then the behavior of observables need to be studied in the limit. 
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This approach was adopted in a recent work on QCD (2+1) by Chakrabarti and 



Harindranath They solved the eigenvalue problem for a few formulations of simi- 

larity renormalization for Hamiltonian matrix elements. To avoid any confusion, we wish 
to emphasize that we use similarity renormalization for particles, and these two frameworks 
can have very different challenges. With this in mind, we note with interest that they find 
convergent eigenvalues, and some structure around a; = 1/2 which they interprete as a sign 
of slow convergence. 

It may turn out that removing the divergences at the level of Hamiltonian is necessary, 
or preferable. There are some compelling reasons for this. Leaving aside the open question 
of whether the limit of bound state solutions with finite 6 converges in the limit 6^0, 
improving our understanding (or lack of thereof) of the vicinity of the zero mode could 
be very useful in constructing more convenient similarity transformations. In present work 
we used as our diagonal part in the similarity transformation the current mass operator. 
Consequently, all formfactors are function of momenta compared to current mass. Even 
though in the bound state equation the kinetic energy is expressed in terms of effective, 
or constituent, A-dependent mass, any expansions in momenta have to be performed in 
reference to current mass. This is certainly not an issue in case of heavy quarks, but for 
light quarks it practically rules out a simple nonrelativistic description of the bound state. 
It is not necessary, though, to use a A-independent diagonal part to drive the similarity 
transformation. It could be devised around a A-dependent, constituent mass, providing we 
have a resolution to the issue of the IR divergences at the Hamiltonian level. 

One possibility is that understanding of the light-cone zero mode is indeed crucial for 
finding the hadronic spectra, and that the requirement that the physical observables are in 
compliance with Lorentz invariance is not sufficient to fully define the IR cutoff theory and 
render it regularization independent. If this is the case, then we are in trouble. 

Another option that comes to mind is more down-to-earth. In a fully covariant theory, 
the physical state surely contains all Fock components. To maintain/restore symmetries 
order by order in coupling, it might be necessary to include in the state at least those Fock 
components that can mix with qq at the given order. To second order the only additional 
Fock component to include would be qqg. It is possible that the infrared divergence in the qq 
cancels with an infrared divergence in qqg. To support this view, recall, that in the standard 
perturbation theory the infrared divergence cancels out between instantaneous exchange and 
one gluon exchange. Even though our framework differs from the usual perturbation theory, 
a similar cancelation might occur. The reason is that the formfactor structure of the mixing 
between qq and qqg is the same as that of the effective qq terms which are plagued by 
the infrared divergence. The main difference is that the mixing depends on a three-body 
wavefunction. Nevertheless, it is conceivable that both the qq effective interaction and 
mixing with the three-body Fock component can be made finite by adding and subtracting 
the same IR counterterm, and it is possible that any ambiguities in connection with the 
choice of a specific IR regulator can be used to improve manifest rotational symmetry in the 
qq sector. 
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